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We  report  gate  bias  dependence  of  the  charge  due  to  piezoelectric  polarization  obtained  by  using  a 
fully  coupled  formulation  based  upon  the  piezoelectric  constitutive  equations  for  stress  and  electric 
displacement.  This  formulation  is  significant  because  it  fully  accounts  for  electromechanical 
coupling  under  the  constraint  of  global  charge  control.  The  coupled  formulation  results  in  lower 
charge  due  to  piezoelectric  polarization  as  compared  to  the  uncoupled  formulation  for  a  given  A1 
mole  fraction.  With  increasing  two  dimensional  electron  gas  concentration,  that  is,  for  gate  biases 
greater  than  threshold,  the  compressive  strain  along  the  c  axis  in  the  barrier  AlGaN  layer  increases 
with  a  concomitant  increase  of  in-plane  stress.  Current  collapse  is  correlated  to  the  increase  in 
source  and  drain  resistances  through  their  dependence  upon  surface  charge.  An  alternate  explanation 
of  current  collapse  using  local  charge  neutrality  is  also  presented.  ©  2006  American  Institute  of 
Physics.  [DOI:  10.1063/1.2203739] 


GaN-based  heterojunction  field  effect  transistors 
(HFETs)  are  currently  being  vigorously  pursued  for  possible 
applications  in  high  temperature  and  high  power  environ¬ 
ment.  Devices  have  been  fabricated  with  superb  power  per¬ 
formance,  notably  the  recently  reported  12  W/mm  at 
10  GFIz  with  a  GaN /  AlGaN/ GaN  high  electron  mobility 
transistor  (HEMT)  without  surface  passivation.1  Further  en¬ 
hancement  and  optimization  of  device  performance  may  be 
possible  with  a  clearer  understanding  of  the  underlying  phys¬ 
ics  of  GaN-based  FIFETs.  In  this  letter,  we  investigate  the 
gate  bias  dependence  of  piezoelectric  polarization  and  its 
possible  implications  on  GaN-based  device  performance. 

The  presence  of  spontaneous  polarization  and  strain  in¬ 
duced  polarization  in  AIGaN/GaN  HFETs  has  been  widely 
recognized.  The  basic  HFET  material  structure  consists  of  a 
thick  GaN  layer  on  which  is  grown  a  pseudomorphic  AlGaN 
layer.  Both  GaN  and  AlGaN  are  grown  in  the  hexagonal 
wurtzite  crystal  structure  with  the  c  axis  normal  to  the  inter¬ 
face.  For  purposes  of  analysis,  the  layers  are  usually  assumed 
to  be  of  infinite  extent  in  the  c  plane,  the  c-plane  lattice 
constant  of  the  pseudomorphic  AlGaN  is  assumed  clamped 
at  the  lattice  constant  of  the  GaN  layer,  and  shear  stresses  are 
assumed  absent.  Because  the  bulk  c-plane  lattice  constant  of 
AIN,  3.112  A,  is  smaller  than  that  of  bulk  GaN,  3.189  A,  the 
clamped  AlGaN  lattice  results  in  tensile  stress  and  strain  in 
the  c  plane  and  compressive  strain  along  the  c  axis  which  is 
assumed  to  be  stress  free.  Following  the  method  used  for 
weakly  piezoelectric  materials,  most  analyses  use  these  as¬ 
sumptions  to  compute  the  c-plane  stresses  and  c-axis  strain 
from  the  tensor  formulation  of  Hooke’s  law,  and  then  com¬ 
pute  the  piezoelectric  polarization  charge  from  the  stress  or 
strain  components  using  the  piezoelectric  moduli.  The  piezo- 
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electric  polarization  is  then  added  to  the  fixed  spontaneous 
polarization  to  perform  further  analysis  of  the  HFET 
operation. With  both  AlGaN  and  GaN  being  rather  strong 
piezoelectric  materials  a  more  complete  approach  would  use 
fully  coupled  electromechanical  equations  to  relate  stress  to 
the  electric  field  as  well  as  to  the  strain  in  these  systems.5  In 
this  letter,  we  apply  the  fully  coupled  equations  to  the 
AIGaN/GaN  HFET.  The  coupled  formulation  results  in  pi¬ 
ezoelectric  polarization  charge  that  is  6.0%  lower  than  the 
simplified  uncoupled  model  for  the  technologically  impor¬ 
tant  AlojGaQjAl/GaN  structure. 

The  coupled  formulation  is  based  on  the  linear  piezo¬ 
electric  constitutive  equations  for  stress  and  electric 
displacement,6 

aij  =  Cijklskl  ~  ekijEk> 

(1) 

Di  =  eijkejk  ~  XijEj  +  P'h 

where  cr/;  is  the  stress  tensor,  Ci;jH  is  the  fourth  rank  elastic 
stiffness  tensor,  skl  is  the  strain  tensor,  etjk  is  the  third  ranked 
piezoelectric  coefficient  tensor,  /c/;-  is  the  second  rank  permit¬ 
tivity  tensor,  I),  is  the  electric  displacement,  Ek  is  the  electric 
field,  and  P\  is  the  spontaneous  polarization.  The  indices  i,  j, 
k,  and  /  run  over  the  Cartesian  coordinates  1,  2,  and  3  (x,  y, 
and  z).  Einstein’s  rule  of  summation  over  repeated  indices  is 
implied.  The  symmetry  of  the  wurtzite  crystal  structure  of 
GaN  and  AlGaN  reduces  the  number  of  independent  elastic 
and  piezoelectric  moduli.  In  the  devices  considered  here,  the 
crystals  are  grown  with  the  c  axis  normal  to  the  surface  in 
the  z  direction,  we  make  the  common  assumption  that  the 
thick  GaN  layer  is  unstrained  and  the  biaxial  strain  of  the 
thin  AlGaN  layer  satisfies  sJC=Sy=(u:Qa^— a AiGaN) AlGaN’ 
where  aGaN  and  a  AlGaN  are  the  c-plane  lattice  constants  of 
each  material.  The  absence  of  stress  along  the  growth  direc- 
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FIG.  1 .  Charge  due  to  piezoelectric  polarization  as  a  function  of  aluminum 
mole  fraction  for  a  pseudomorphic  AlGaN  layer  on  GaN. 


tion  (z  direction),  in  the  barrier  AlGaN  layer,  allows  us  to 
express  the  strain  along  the  growth  direction  as 


■2^e,+  ^£A1GaN 

C33  C33  z 


(2) 


where  FA1GaN  is  the  ^-directed  electric  field  in  the  barrier 
AlGaN  layer  and  we  have  expressed  the  elastic  and  piezo¬ 
electric  moduli  in  matrix  notation.7 

The  areal  charge  concentration  due  to  piezoelectric  po¬ 
larization  in  AlGaN  is  expressed  as  follows: 

nAlGaN  0  (  ^13  \  AlGaN  e33 

^PE  =  2sd  e3!  -  —  e33 1  +EZ  — .  (3) 

\  c33  /  c33 


It  should  be  noted  that  in  uncoupled  analyses,  the  second 
term  on  the  right  hand  side  of  Eq.  (3)  is  absent. 

Applying  the  continuity  of  perpendicular  components  of 
the  electric  displacement  vector  at  the  barrier  AlGaN  /GaN 
heterointerface,  we  express  the  field  in  AlGaN  as  follows: 


.AlGaN 


_ 1 _ 

KA1GaN  +  (4/c33) 


cr2D  +  A  Ps 


e3r 


Cp 

C 


e33 


33 


(4) 


where  /cAIGaN  is  the  permittivity  of  the  barrier  AlGaN  layer 
and  AE'  =  EGaN_^AiGaN  is  the  difference  between  the  spon¬ 
taneous  polarization  between  GaN  and  AlGaN.  The  deriva¬ 
tion  of  the  Eq.  (4)  is  based  upon  the  assumption  that  the 
unintentional  doping  concentration  of  GaN  is  negligible. 
This  assumption  allows  us  to  relate  the  electric  displacement 
in  GaN,  at  the  AlGaN  /GaN  heterointerface,  to  the  two 
dimensional  electron  gas  (2DEG)  charge  density,  <x2D,  as 
^GaNgGaN  s  Substituting  (4)  in  (3)  gives  the  piezoelec¬ 
tric  component  of  total  polarization  in  AlGaN  as 


<GaN  =  2ex(e3l  -  ^Je33|(l  -  «)  +  «(cr2D  +  AEO,  (5) 

where  a=(e33/C33)/(/cA1GaN  +  e33/C33)  is  a  measure  of  the 
electromechanical  coupling.  Note  that  in  the  absence  of  elec¬ 
tromechanical  coupling  a= 0  and  we  revert  back  to  the  ex¬ 
pression  obtained  by  applying  Hooke’s  law  and  piezoelectric 
polarization  in  an  uncoupled  way. 

Figure  1  shows  plots  of  the  room  temperature  areal 
charge  concentration  due  to  piezoelectric  polarization  with 
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FIG.  2.  Piezoelectric  charge  component  for  a  pseudomorphic  Al^Ga^N 
layer  on  GaN  as  a  function  of  the  2DEG  concentration  with  aluminum  mole 
fraction,  x,  as  a  parameter.  The  solid  lines  represent  the  variation  of  piezo¬ 
electric  charge  density  using  the  modified  formulation  and  dotted  lines  show 
the  uncoupled  charge  density.  Dashed  lines  show  the  corresponding  electric 
field  in  fully  coupled  case. 


and  without  electromechanical  coupling  for  the  case  of  neg¬ 
ligible  2DEG  concentration  [cr2D=0  in  Eq.  (5)],  which  would 
occur  for  a  HFET  gate  biased  below  threshold.  For  compu¬ 
tations  we  used  the  material  constants  reported  in  Ref.  4.  The 
absence  of  electromechanical  coupling  results  in  an  over  es¬ 
timation  of  the  charge  due  to  piezoelectric  effect.  From  Eq. 
(5),  we  see  that  there  are  two  factors  contributing  to  the 
lower  areal  charge.  First,  the  polarization  due  to  the  strain  is 
reduced  by  a  factor  of  ( 1  -  a) .  When  electromechanical  cou¬ 
pling  is  included,  the  amount  of  stress  in  the  c  plane  of  the 
AlGaN  layer  is  lower,  resulting  in  lower  piezoelectric  polar¬ 
ization.  Second,  the  coupling  of  the  electric  field  due  to 
spontaneous  polarization  has  the  effect  of  further  relaxing  the 
tensile  stress  of  the  AlGaN  layer,  again  lowering  the  piezo¬ 
electric  polarization. 

The  application  of  gate  bias  greater  than  threshold  in¬ 
duces  a  2DEG  concentration,  cr2D,  in  the  channel  of  the 
HFET.  From  Eq.  (5)  it  is  clear  that  the  change  in  cr2D  is 
accompanied  by  a  change  in  the  piezoelectric  charge  in  Al¬ 
GaN  layer.  Figure  2  shows  the  areal  charge  concentration 
due  to  the  piezoelectric  effect  as  a  function  of  2DEG  con¬ 
centration  with  A1  mole  fraction  as  a  parameter  for  both  the 
coupled  and  the  uncoupled  analyses.  When  electromechani¬ 
cal  coupling  is  neglected,  the  polarization  charge  is  indepen¬ 
dent  of  cr2D.  When  coupling  is  included,  the  areal  charge 
increases  with  increasing  2DEG  concentration  for  a  given  A1 
mole  fraction,  moving  toward  the  uncoupled  value. 

With  the  increasing  2DEG  concentration  corresponding 
to  the  applied  gate  bias  becoming  less  negative,  the  compres¬ 
sive  strain  along  the  c  axis  in  the  barrier  AlGaN  layer 
increases.5  Figure  3  shows  the  c-axis  strain  of  an  Al0  3Gao  7N 
barrier  layer  as  a  function  of  2DEG  concentration  with  and 
without  the  coupled  electromechanical  formulation.  In  the 
uncoupled  formulation,  the  lattice  is  compressively  strained 
by  a  constant  value  of  0.386%.  However,  with  the  incorpo¬ 
ration  of  the  electromechanical  coupling,  the  magnitude  of 
the  strain  in  the  barrier  layer  increases  from  0.318%  for  a 
2DEG  concentration  of  1  X  1012  cm-2  to  0.358%  for  a  2DEG 
concentration  of  1  X  1013  cm'2. 
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FIG.  3.  Coupled  strain  and  uncoupled  strain  in  the  c-axis  direction  of  a 
pseudomorphic  Al^Ga^N  layer  on  GaN  as  a  function  of  2DEG  concentra¬ 
tion  with  aluminum  mole  fraction,  x,  as  a  parameter. 


To  account  for  the  fact  that  the  spontaneous  and  piezo¬ 
electric  polarizations  produce  dipoles  with  inherently  bal¬ 
anced  charges,  we  have  implicitly  assumed  the  presence  of 
surface  and  substrate  charges  that  are  responsible  for  the 
2DEG  concentration.  Using  global  charge  neutrality,  the 
2DEG  concentration  can  be  expressed  as  the  difference  of 
the  surface  and  substrate  charge  concentrations  allowing  us 
to  recast  Eq.  (5)  in  the  following  form: 

^reGaN  =  2s_t|e31  -  ^e33j(l  -  a)  +  «(o-surface  +  AP!), 

(6) 

where  we  have  neglected  the  contribution  of  the  substrate 
charge.  Global  charge  neutrality  requires  that  any  modifica¬ 
tion  in  the  surface  charge  be  accompanied  by  a  matching 
change  in  the  2DEG  concentration.  The  role  of  polarization 
is  in  the  modification  of  the  electric  field  in  the  AlGaN  bar¬ 
rier  layer,  as  shown  in  Fig.  2.  With  an  increase  in  the  positive 
surface  charge,  PpgGaN  increases  with  a  concomitant  decrease 
in  the  electric  field  as  well  as  a  corresponding  change  in  the 
surface  potential;  this  in  turn  must  be  consistent  with  the 
surface  charge. 

The  2DEG  and  the  associated  parasitic  gate-source  and 
gate-drain  resistances  are  thus  controlled  by  the  surface 
charge.  For  example,  applying  surface  passivation  to  modify 
the  surface  charge  in  the  gate-source  and  gate-drain  regions 


will  alter  the  source  and  drain  resistances.  Modification  of 
the  surface  charge  in  the  gate-drain  and  gate-source  regions 
may  also  occur  when  a  negative  gate  bias  causes  electrons  to 
be  injected  from  the  drain  and  source  ends  of  the  gate  result¬ 
ing  in  a  decrease  in  the  net  positive  surface  charge  in  those 

g 

regions.  This  lowers  the  induced  charge  at  the  AlGaN  /  GaN 
heterointerface,  increasing  the  parasitic  drain  and  source  re¬ 
sistances.  The  increase  in  parasitic  resistances  due  to  surface 
charge  may  be  a  mechanism  of  current  collapse.  As  an  alter¬ 
native  to  global  charge  neutrality,  local  charge  neutrality  can 
be  used  to  equate  the  2DEG  concentration  to  the  difference 
in  polarization  charges  at  the  AlGaN  /GaN  heterointerface. 
In  this  interpretation,  a  reduction  in  PpgGaN  [Eq.  (6)]  can 
attributed  to  the  relaxation  of  in-plane  stress  in  AlGaN  layer, 
in  agreement  with  the  observations  made  by  Simin  et  al?  It 
should  be  noted  that  the  dynamics  of  current  collapse  will 
depend  upon  the  type  of  traps  and  their  associated  trapping 
and  detrapping  time  constants. 

In  conclusion,  we  have  presented  a  coupled  electrome¬ 
chanical  formulation  for  the  self-consistent  determination  of 
piezoelectric  charge  in  AlGaN/GaN  heterointerfaces.  The 
uncoupled  formulation  overestimates  the  charge  associated 
with  piezoelectric  polarization  by  overestimating  the  stress  in 
the  c  plane  and  neglecting  the  field  due  to  spontaneous  po¬ 
larization.  Any  surface  modification  caused  by  passivation  or 
electron  injection,  for  example,  will  change  the  carrier  con¬ 
centration  in  the  neighborhood  of  AlGaN/GaN,  thereby 
changing  the  source  and  drain  resistances  and  influencing 
current  collapse. 
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